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Abstract 
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I. INTRODUCTION 



The simplest modified Morse potential model suggested by Deng and Fan [1] and related 
to the Manning- Rosen potential [2] (also called Eckart potential by some authors [3]) is 
anharmonic potential defined by 



where r e (0, oo), and the three positive parameters and a denote the dissociation en- 
ergy, the equilibrium inter-nuclear distance and the range of the potential well, respectively. 
The above potential is used to describe diatomic molecular energy spectra and electromag- 
netic transitions and is the true internuclear potential in diatomic molecules with the same 
behaviour for r — )■ [4]. The above potential was called a generalized Morse potential (GMP) 
model and illustrated in Figure 1 for various values of potential parameters. As stated in 
Ref . [4] , the Morse potential and the GMP are very close to each other for large values of Ve 
{a — 1) in the regions r ^ re and r > but are very different at r ~ 0. Further, if the two 
potentials are deep (D 3> 1), they could be well approximated by a harmonic oscillator in the 
region r ~ rg (see Fig.l in [4]). To describe the vibrational spectra of diatomic molecules the 
potential curve V{r) is approximated by a sum of three Morse functions [5]. The approach 
has the advantage of being more flexible than the simple Morse potential while preserving 
the correct asymptotical V{oo) — const [5]. The potential model (1) is a special case of the 
five-parameter exponential-type potential model [6,7]. The exact solvability of the s-wave 
(Z = case) bound state energy eigenvalues and eigenfunctions of the GMP is due to the fact 
that it belongs to the class of the Eckart potential, a member of the hypergeometric Natan- 
zon potentials which can be solved algebraically by means of S0{2, 2) symmetry algebra [4] 
and SO {2, 1) algebra [8]. The X-H stretching motion in small molecules has been treated by 
the potential model (1) [9]. Moreover, the approximated Z-wave bound state solutions [10] 
of the Schrodinger equation has been solved using the conventional approximation scheme 
suggested by Greene and Aldrich [11] to deal with the centrifugal barrier term l{l + l)/r^, 
singularity at r = 0. The exact analytic expressions for matrix elements of positive integral 
powers of the coordinate [12] and the quasi-one-dimensional system of DNA [13] have also 
been studied. 

To study the relativistic effects and corrections of the molecular Morse potential [4] , the 
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Dirac equation has been solved for attractive scalar S{r) and repulsive vector V{r) Morse 
potentials under pseudospin symmetry in the nuclear theory using the Pekeris approximation 
[14]. Recently, the approximate bound state solutions of the pseudospin and spin symmetric 
Dirac equation with the GMP has been calculated using an improved approximation scheme 
to deal with the centrifugal (pseudo-centrifugal) term and by employing the basic concept 
of the supersymmetric shape invariance formalism (cf. [15] and the references therein). 

Many authors have investigated approximately the solution of the Dirac equation with 
a few potential models such as the Hulthen potential [16], the Hulthen potential including 
Coulomb-like tensor potential [17], the generalized Woods-Saxon potential [18,19], the Eckart 
potential [20], the Morse potential [14,21], the Poschl-Tcller potential [22], the Manning- 
Rosen potential [2,23], the hyperbolic potential [24], the Rosen-Morse potential [25], the 
pseudoharmonic potential [26] , and the Kratzer potential connected with an angle-dependent 
potential [27], etc within the framework of various methods. 

The diatomic molecular model consisting of nuclei having masses mi and m2, the reduced 
mass is defined /i — 777.1777.2/ (mi -|- m2) can be included to the spin symmetry and the 
pseudospin symmetry concepts [28]. Ginocchio [29-31] showed that the spin symmetry occurs 
when the difference potential between the vector potential V{r) and scalar potential S{r) is 
a constant (i.e., A(r) = V{r) — S{r) = Cs — constant) and the pseudospin symmetry occurs 
when the sum potential of the vector potential V{r) and scalar potential S{r) is a constant 
(i.e., E(r) = V{r) —S{r) — Cps — constant). The spin symmetry concept [32] is particularly 
relevant for mesons [33]. The pseudospin symmetry concept [34,35] in nuclear physics refers 
to the quasi-degeneracy of single-nucleon doublets and can be characterized with the non- 
relativistic quantum numbers (n, /,j = / + 1/2) and (n — 1, / + 2, j = / -|- 3/2), where 77, / 
and j are the single-nucleon radial, orbital and total angular momentum quantum numbers, 
respectively. Alhaidari et al. [36] investigated in detail the physical interpretation on the 
three-dimensional Dirac equation in the cases of spin symmetry limitation A(r) = and 
pseudospin symmetry limitation E(r) = 0. In real nuclei, E(r) 7^ constant and pseudospin 
symmetry is only an approximation. The quality of the pseudospin symmetry approximation 
depends on the competition between the pseudo-centrifugal potential and the pseudospin 
orbital potential [37]. 

Recently, in the framework of the spin symmetry S{r) ~ V{r) and pseudospin symmetry 
S{r) ~ — y(r), the bound state energy eigenvalues and associated upper- and lower-spinor 
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wave functions are investigated by means of the Nikiforov-Uvarov (NU) method [38]. We 
have approximately solved the Dirac equation for the Rosen-Morse potential [25] with spin 
and pseudospin symmetry for any k state and found the eigenvalue equation and corre- 
sponding two-component spinors within the framework of an approximation to the term 
proportional to We have also solved the (3-1-1) dimensional Dirac equation for a par- 

ticle trapped in the spherically symmetric generalized WS potential under the conditions of 
exact spin and pseudospin symmetry combined with approximation for the spin-orbit cen- 
trifugal (pseudo-centrifugal) term, and calculated the two-component spinor wave functions 
and the energy eigenvalues for any arbitrary spin-orbit k bound states [18]. 

In principle, the solution of the Schrodinger equation for the GMP model can be used to 
describe the motion of the nucleons in the mean field produced by the interactions between 
nuclei. However, the Dirac equation successfully merges quantum mechanics with special 
relativity and is considered to be the natural transition to quantum field theory. It provides 
a natural description of the electron spin, predicts the existence of antimatter and is able 
to reproduce accurately the spectrum of the hydrogen atom. It also predicts some peculiar 
phenomena such as Klein's paradox and unexpected quivering motion of free relativistic 
quantum particle which are key examples for understanding relativistic quantum effects, 
but are difficult to observe in real particles. In order to to understand the origin of spin 
and pseudospin symmetry, we need to take into consideration the motion of the nucleons 
in a relativistic mean field and consider the Dirac equation [31]. We attempt to study the 
approximate solution of the Dirac equation for the GMP with non-zero spin-orbit quantum 
numbers k by employing an improved approximation scheme (see e.g., [17,39,40] and the 
references therein) to deal with the spin-orbit centrifugal, k{k, + l)/r^ (pseudo-centrifugal, 
k(k — l)/r^) barrier term. The inter-relation between the GMP with the Oscillator and 
Morse potentials investigated in Ref. [4] and the recent relativistic treatments of the Morse 
and Oscillator potentials in [14,26] are some motivations for the present study. The simple 
transformation of the Dirac equation into the Schrodinger-like equation and the success in 
studying the approximated bound state solutions of the Dirac equation with various potential 
models [14,17,18,25,26] also make the solution possible. 

Wc tend to show that the new scheme of parametric generalization of the NU method 
[41] given in Appendix A is a powerful tool for solving a second order differential equation 
by turning it into a hypergeometric type equation [38]. The advantage of employing the 
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NU method, in the present work, is that it can be used to find the bound state energy 
spectra and the corresponding spinor wave functions under the condition of spin symmetry 
and pseudospin symmetry concept for any k state in a very simple way. 

This paper is organized as follows. In section 2, we investigate the bound state energy 
equation and the corresponding two-component spinor wave functions under the condition 
of spin symmetry and pseudospin symmetry concept for the GMP model by employing a 
parametric generalization of the NU method. In section 3, we study some special cases like 
the s(S)-wave cases k — ±1 {I — I — 0), the non-relativistic limit and the a — >■ (Kratzer 
potential). In section 4, we present some numerical results to the non-relativistic and rela- 
tivistic numerical energy levels for GMP and Kratzer models. The relevant conclusions are 
given in section 5. 




II. DIRAC BOUND STATE SOLUTIONS 

The Dirac equation for fermionic massive spin- 1/2 particles moving in an attractive scalar 
potential S{r) and a repulsive vector potential V{r) is given by [42] 

[ca • p + /3 (Mc^ + S{r)) + V{r) - E] V^„,(r) = 0, V'^Jr) = V^(r, 9, cj>), (2) 

where E is the relativistic energy of the system, m is the mass of a particle, p = —iTiV is 
the momentum operator, and a and (3 are 4x4 Dirac matrices, i.e., 

I0\ /01\ /O-A /lO. 

(3) 

where I denotes the 2x2 identity matrix and (Tj are the three-vector Pauli spin matrices. 
For a spherical symmetrical nuclei, the total angular momentum operator of the nuclei J and 
spin-orbit matrix operator K = — /3 (cr • L + I) commute with the Dirac Hamiltonian, where 
L is the orbital angular momentum operator. The spinor wave functions can be classified 
according to the radial quantum number n and the spin-orbit quantum number k and can 
be written using the Pauli-Dirac representation in the following forms: 

where the upper- and lower-spinor components Fn^{r) and Gn^{r) are the real square- integral 
radial wave functions, and the spin spherical harmonic functions 



coupled to the total angular momentum j and it's projection m on the z axis and k{k + 1) — 
l{l + 1) and k{k — 1) — l{l + 1). The quantum number k is related to the quantum numbers 
for spin symmetry / and pseudospin symmetry I as 

+ = + (si/2,P3/2, ete.), j = ^ + aligned spin < 0) , 
+/ = + (j + \) , (pi/2, d3/2, etc.), j unaligned spin > 0) , 

and the quasi-degenerate doublet structure can be expressed in terms of a pseudospin angular 
momentum 's — 1/2 and pseudo-orbital angular momentum / which is defined as 



K — 



-1 = - {j + ^) , {si/2, P3/2, etc.), 3=1- 1/2, aligned spin (k < 0) , 
+ iT+ ^ = + + i) , (c?3/2, /5/2, etc.), j = T+ 1/2, unaligned spin > 0) , 

(6) 

where n = ±1,±2, ••• . For example, (lsi/2, 0^3/2) and (2p3/2, 1/5/2) can be considered as 
pseudospin doublets. 

Upon direct substitution of Eq. (4) into Eq. (2), we can obtain two radial coupled Dirac 
equations for the two spinor components as follows: 

{'t + 7) ^"'^^^^ " t^'^' + ^"'^ " ^^^^ 

{'i ~ 7) ^"'^^''^ ^ ~ ^"-^ + ^^""^^ ^^^^ 

where A(r) = V{r) — S{r) and E(r) = ^(r) -|- 5'(r) are the difference and sum potentials, 
respectively. 

Under the spin symmetry ( i.e., A(r) = C5 = constant), one can eliminate Gnnif) in 
Eq. (7a), with the aid of Eq. (7b), to obtain a second-order differential equation for the 
upper-spinor component as follows: 

(f k{k+1) 1 



Fn.(r) 



= ^ {El^ - M'c' + (Mc^ - K.)) i^n«(r), (8) 
and the lower-spinor component is obtained from Eq. (7a) as 

where 7^ —M(?, only real positive energy states exist when = (exact spin symme- 
try). On the other hand, under the pseudospin symmetry ( i.e., E(r) = Cpg = constant), 
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one can eliminate Fnnir) in Eq. (7b), with the aid of Eq. (7a), to obtain a second-order 
differential equation for the lower-spinor component as follows: 

1 



(Mc2-£;„« + Cp,)A(r) 



1 



and the upper-spinor component Fni^{r) is obtained from Eq. (7b) as 

1 



ps 



(10) 



(11) 



where Enn 7^ Mc^, only real negative energy states exist when Cps = (exact pscudospin 
symmetry). It is worthy to note that the reality and finiteness of our solutions demand that 
the upper and lower radial components should satisfy the essential boundary conditions: 
FnM = ^^^(O) = and F„«(oo) = ^^^(oo) = 0. 



A. Spin symmetry solutions of the GMP model 



At first, we investigate the spin symmetry by taking the S(r) = 2V{r) — )■ VcMpir) as 
mentioned in Ref. [43] enables one to reduce the resulting relativistic solutions into their 
non-relativistic limit under appropriate transformations. From Eq. (8), wc can sec that the 
energy eigenvalues, EnK, depend only on n and I, i.e., E^k — E{n,l{l + 1)). For I ^ 0, the 
states with j — I ± 1/2 are degenerate. This is a SU{2) spin symmetry. Following Refs. 
[22-25], we impose the GMP [1,4] as the E(r), i.e., 

^ir)-D(l--±-)\ (12) 



e"' — 1 / 

leads us to obtain a Schrodinger-like equation in the form: 

d'^ k(k+1) oof-. &e-"^ ^ ^ 
dr^ V 1 - e-" 



I 2 2 
+ a 



Fnn{r) = 



(13) 



1 



{Mc^ + E^,-C,) D, 



1 



[EL-^'c^ + C.(Mc2-S„.)], (14) 



where k (k -|- 1) = Z (Z -|- 1) , k = Z for k < and k = — (Z -|- 1) for k > 0. The exact solution of 
the above equation is possible only for the s-wave case {k — —1) due to the centrifugal term 
K (k -I- 1) jr^. However, if k is not too large, the case of the vibrations of small amplitude 



about the minimum, we attempt to use the following improved new approximation scheme 
to deal with the centrifugal (pseudo centrifugal) term, near the minimum point r — re {i.e., 
x^O), (cf. Refs. [17,40,44-48]): 

n2 



a 



do + 



:i-e- 



a 



do + 



(ar 



1 {ary (ar)^ (ar)^ 
12 ^ 240 ~ 6048 ^ 172800 



+ O {{arf) 
(15) 



When ar <C 1, the value of the dimensionlcss constant do = 1/12 has simply determined by 
the above series expansion and a takes the unit of reciprocal of r. The present approximation 
was shown to be more powerful than the usual approximation [40]. Obviously, the above 
approximation to the centrifugal (pseudo-centrifugal) term turns to 1/r^ when the parameter 
a goes to zero as 



lim 



aM rfo + \ o I = (16) 

which shows that the usual approximation is the limit of our approximation (cf. e.g., [17,40] 
and the references therein). 

We introduce the following new dimensionlcss parameter, z(r) — e""'' G [0, 1], which 
maintains the finiteness of the transformed wave functions on the boundary conditions. 
Thus, substituting Eq. (15) into Eq. (13), we obtain the following Schrodinger-like equation 
satisfying F^nir), 

d^ (l-z) d \ ^ . , 



+ 



z\\-zf 
where 



dz^ z{\ — z) dz ^ 

{- [(2 + h) hv\ + + [2 {hv\ + 4J _ ^ + 1)] ^ _ 4J F„«(^) = 0, 

(17) 



el^ = I/? - a;? + «;(« + 1) do, (18) 

and Fnnir) = FnK{z). In order to clarify the parametric generalization of the NU method 
[25], let us take the following general form of a Schrodinger-like equation written for any 
potential. 



+ 



dz"^ a{z) dz (^"^{z, 



satisfying the wave functions 



where 



V'n(^) = (t>{z)yn{z), 
t{z) = Ci - C2Z, 



(19) 

(20) 
(21) 



and 

a{z) = ^ (1 - c^z) and a{z) = -Az^ + Bz-C, (22) 

are two polynomials at most of first- and second-degree, respectively. Furthermore, when 
Eq. (17) is compared with its counterpart Eq. (19), we can obtain the specific values for 
the constants Cj {i — 1,2, 3) along with A, B and C. Now, following the NU method [38] 
and making the substitution of Eqs. (21) and (22), we can obtain general forms for the 
polynomials 7r{z) and t{z), the root of the parameter k, the eigenvalues equation and the 
wave functions (l>{z) and yn{z) expressed in terms of the constants q (i = 4, 5, • • • , 13) as 
shown in Appendix A (cf. Refs. [18,21,25,40]). Therefore, the task of computing the energy 
eigenvalues and the corresponding wave functions of Eq. (17) within the framework of 
the parametric generahzation of the NU method is relatively easy and straightforward. It is 
explained in the following steps. Firstly, we need to find the specific values for the parametric 
constants Cj (i = 4, 5, • • • , 13) by means of the relation Al of Appendix A. The values of 
all these constants q (i = 1, 2, ■ ■ ■ , 13) together with A, B and C are therefore displayed in 
Table 1 for the GMP model. Secondly, using the relations (A2-A5), the analytic forms of 
the essential polynomials t:{z) and t{z) along with the root A;, required by the method [41], 
can also be found as 



z I 
= -2 - 2 



{j{l + 2Kf + Ah^ul + 2e^^ - 2£, 



(23) 



k = 26i/2 _ ^ + 1) _ £^^^(1 + 2kY + 46V2, (24) 

t{z) = 1 + 2£„, - 2 (^1 + + ^^(l + 2«;)' + 46V?^ z, (25) 

where r'{z) < must be satisfied in order to obtain a physical solution according to the NU 
method [38] . Thirdly, we need to calculate the energy eigenvalues by means of the eigenvalue 
equation, relation A6, and obtain 

_ {2 + h)hul 2n + l + -^(l + 2/t)V46V^ 

£n« — I . ■ (26) 

2n + l + -^(l + 2K)^ + 46V2 ^ 

Finally, with the aid of Eqs. (14) and (18), Eq. (26) can be also reduced to the energy 

equation for the GMP with the spin symmetry concept for any spin-orbit quantum number 
K — ±1, ±2, • • • values in the Dirac theory, 

(Mc^ + - Cs) {Mc^ + D- K«) + h^cWK {k + 1) do 



2(n + 5i) 



, n = 0,l,2,--. , 



(27) 



where 



5i = ^ (^1 + ^(l + 2«)' + 46V?^ > 1. (28) 

In what follows, in order to establish the wave functions Fn^ir) of Eq. (8), the relations 
(A7-A10) are used. Firstly, we find the first part of the wave functions yields 



(P{Z)=Z'-{1-ZY\ Snn>0, ^1 > 0. 

Secondly, we calculate the weight function as 

and this, in turn, generates the second part of the wave functions. 



(29) 



(30) 



l/„(^)~z-^^"«(l-z)-(^^-i) — 



n+25i-l 



p^2e„.,25i-l)^^_2^)^ (31) 



where Pji^'^\l — 2z) is the orthogonal Jacobi polynomials [49,50]. Finally, the upper spinor 
component Fn^{z) for arbitrary k can be obtained by means of Eq. (20) as 



_ r{n + 2en. + l) r 



Fnnir) = A/'n«e-""'^"^(l - e-"^)5ip^2£„.'25i-i)^^ _ 2e-ar-) 

1 - e-"0^i (-n, n + 2£„« + 25i; 1 + 2£„,; e""'') , (32) 



r(2£„, + l)n\ 

where the normalization constants Hnn are calculated in Appendix B. 
Let us recall the derivative relation of the hypergeometric function, 

ah 



d 

dz 



iFi {a;b;c;z) 



,Fi {a + l;b + l;c+l;z), 



that is used to calculate the corresponding lower-component Gnnir) in Eq. (9) as 

J^UK 



Mc^ + - V 1 - e 



— - asnn + - F,,^{r) 



,*/- na (n + 2£„^ + 2(5i) c>r^^l / -ar\^n.+l 



) (1 + 2£„k) 



X 



'(Mc2 + 

(1 - n; n + 2 (e,, + ^0 + 1; 2 (1 + £,,) ; e""^) . 



(33) 



We would like to note that the hypergeometric series 
2F1 (1 — n; n + 2 (£„K + <^i) + 1; 2 (1 + Sns) ; e""'') terminates for n = and thus does 
not diverge for all values of real parameters and 5i. 
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In the presence of exact spin symmetry (Cg — 0), Enn 7^ —Mc^ (only positive energy 
states do exist). In the exact spin symmetry [Cg — 0) with {k — 1,k — —2), the upper- 
and lower- spinor wave functions for the ground (Opi/2, 0^3/2) and first excited (lpi/2, IP3/2) 
degenerate eigenstates are being illustrated in Figure 2a and 2b, respectively. A glance at 
Figure 2 reveals that there is only one curve (dashed line) for the radial wavefunctions of 
the upper components for both states in the doublet. However, there are two curves (solid 
lines) for the radial wavefunctions of the lower components. The following values of the 
parameters M = 1.0 fm~^, D — lb fm~^, a = 0.1 fmT^ and Ve — 0.4 /m, -Boa=i — 
-Eo,«=-2 = 5.5791076 fm~^ and E^^^^i = -Ei,«=_2 = 8.1823677 fm'^ have been used. 

Let us now study the special case when the parameter a ^ in Eq. (1). The GMP 
potential can be easily reduced to the well-known Kratzcr molecular potential, 

lim y(r) = L> f ^^^^ . (34) 



which has been studied extensively by using different methods as the function analysis [51], 
the NU [52,53] and the exact quantization rule (EQR) [54]. To avoid the repetition, following 
Appendix A, we write down the essential polynomials: 

7r(r) = i[l + 7-2e„,r], (35) 
2r 

A; = (Mc^ + -Cs)D- 7e„„ (36) 

it c 

r(r) = l + 7-2e,,r, (37) 

where 



4r^ 

7 = -y/ (1 + 2kY + ^ (Mc2 + K« - Cs) D, (38) 



e„. = ^ V(Mc2 - + D) {Mc^ + - C,). (39) 

We further obtain the following two expressions which are relevant in the construction of 
the energy equation (relation A6) [38] 

2r 

A„ = 2ne„« and A = -^-^ (Mc' + K« - C,) L> - (1 + 7)e„«, (40) 

ft c 



and substituting = A, we can obtain the following spin symmetric energy equation, 

(Mc^ - + D) ^ ^^'iMc^ + E..-Cg) ^^^^^ 

[Nn + ^/Nl + qDiMd^ + Er^.-CsYj 
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where Nn — 2n + 1, N^, — 2k + 1 and q — {2re/Ticf . The above equation for energies looks 
hke a quartic equation of the form: 

a^E^^ + a^El^ + 02^;^ + «i^n« + oo = (41b) 

with coefficients 

04 = q^D''- as = 2qD - - qDCs] ; 

a2 = {Nl - Nlf + 2qD{D + M + C,)(iV^ - N^) + AqD^Nl + q^D^ (C^ ^ ^MCs - 2M^) ; 

ai = 2q^D'^MCs{M - C,) - 2{D + M)(7V^ _ ^2^2 

^-2qD (M^ - 2MC, - DCs) {K - N^) - AqD^ {D + C,) N^; 



ao = \qDM^ - {D + M){Nl ~ N^)] ' - AqD^ (M - C,) {D + M)N^ 
+q^D^M^Cs {Cs - 2M) + 2qDMCs{D + M){NI - NI)] 

where we have set c = 1. For a given value of n and n (or /), the above quartic equation, Eq. 
(41b), provides four distinct positive and negative real (real and complex) energy spectra 
related with or E~^, respectively. One of the distinct solutions is only vahd to obtain the 
positive-energy bound states in the limit of the spin symmetry. Therefore, the procedures for 
calculating the four distinct energies; namely. En}, En}, En} and En} are given in Appendix 
B. 

Furthermore, following [18], the normalized upper- and lower-spinor wave functions can 
be calculated as 



and 



X 



Mc2 + En^ - Cs ' ' \l {n + K + l)T{n + 2K + 2) 
(^^^^ - Fn.{r) - 2en.r^^'e-^-^L^:^^'\2en.r) 



^ = ^(7-1), (43) 
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respectively, where L^i\x) are associated Laguerre polynomials. The simplest exact spin 
solution, representing the ground state and first excited state, are 

2 



1 



(2eoK 



(2eo«) 



Ko+l 




3+^(l+2«)^+^(Mc2+Bi« 

^i^i+ig-ei«r (_2ei„r + 2X1 + 2) , 



^iiTi+lg-eiKf 



respectively. Finally, we would like to note that eo^ and ei^ can be calculated via Eq. (39) 
whereas Kq and Ki via Eq. (43) along with Eq. (38) when Cg — 0. 



B. Pseudospin symmetry solutions of the GMP model 



From Eq. (10), we can sec that the energy eigenvalues, Enn, depend only on n and /, i.e., 
EnK = E{n, l{l + 1)). For / 7^ 0, the states with j = / it 1/2 are degenerate. This is a SU{2) 
pseudospin symmetry. Following Refs. [22-25], we impose the GMP [1] as the A(r), i.e., 

AW-^(l-i^)'' (44) 
leads us to obtain a Schrodinger-like equation in the form: 



(P k{k-1) ^ ^2^2 ^ 
dr"^ j,2 2 



fee- 



where 



1 



1-e- 



1 



I 2 2 
+ a UJ2 



Gnn{r) = 



(45) 



[£;^^-MV-Cp,(Mc2 + (46) 



where k{k — 1) — l{l + 1). We follow the same procedures in the previous subsection to 



obtain a Dirac equation satisfying Gn^ir), 

f (P (l-z) d 



- 



+ 



ydz"^ z(l — z) dz 
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+ ^2(1^_^)2 (- - (2 + &) &^2] + [2 - &^2) " (^^ " 1)] ^ "^L} Gnn{z) = 0, 

(47) 

where 

= -u;^ - + (/^ - 1) do, (48) 

where Gnnii^) = Gnniz)- To avoid repetition in the solution of Eq. (45), a careful inspection 
for the relationship between the present set of parameters (£„«, 2^2) ^"^^ the previous one 
{SriKi i^i) tells us that the negative energy solution for pseudospin symmetry, where S{r) ~ 
—V{r), can be obtained directly from the spin symmetric solution by using the following 
parameter mapping [18,55]: 

Funir) ^ Gnnir), K ^ K-1, V{r) -V{r) (i.e., D -D), E^^ ^ and Cs -Cps- 

(49) 

Following the previous procedures, the constants in the case of pseudospin symmetry concept 
are displayed in Table 2. Applying the above transformations to Eqs. (27), (32) and (33) 
lead to the following pseudospin symmetric energy equation, 

- Mc^ - Cps) {D - Mc" - K«) + n''c'a\ (k - 1) do 

^,.,„.(|i^^^M)\„.o,.v.., 

with 

5,^^(l + ^{l-2Kf-4b^ui^ >1. 

Note that Eq. (50) can be also expressed in the form of quartic equation (cf. Eq. (41b)) 
since the two parameters 1^2 and 62 contain the energy eigenvalues EnK- The procedures 
of this analytic solution is so similar to the one presented in Appendix B with the changes 
D — > — D, Eni^ — > —EnK and Cg — > —Cps- Furthermore, the lower-component wave functions: 

= ■A4. ^^^"'^"^^^" e-^-"^(l - e'^^y^ (-n, n + 2£„, + 25^; 1 + 2?„«; e""'') , (52) 

77/. 

where 

(2 + b) hul 2n+l + -^(l-2/.)^-4&V| 



(51) 
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The upper-component -F„„(r) can be calculated from Eq. (11) as follows 



K 



MC^ - Enn + Cps 

na (n + 2e„K + 2^2 



{Md'-Enn + Cps){l + 2e, 
X 2F1 (l-n;n + 2(?„„ + 52) + l;2(l + £„«);e-"''). (54) 

Hence, in the exact pseudospin symmetry where Cps = 0, Enn 7^ Mc^ (only negative energy 
states exist). 

On the other hand, the pseudospin solutions of the Dirac equation for the Kratzer po- 
tential can be obtained from the spin symmetry case by applying transformation map given 
by Eq. (49) as 

(55) 



7= W(1-2k)^ 



ADrl 



(Mc''-E,,, + C,,). 



en. - J-\/-D (Mc-^ - + Cp,) - {El^ - M-^c^ - C^, (Mc^ + E^,)). 



(56) 



Therefore, the eigenvalue equation is 



El^ - M^c^ - {Mc' + E„,) = -D {Mc^ - E^, + C^,) 



1 fDre(Mc''-En. + Cps) 



h^c^ V n + K + 1 
and the normalized lower- and upper-spinor wave functions are given by 



(57) 



Gnnir) = (2e„«) 



K+l 



\^ + X + 1) r (^n + 2is: + 2) 



{lEnnT), (58) 



Fr,Ar) 



X 



Mc2 - + Cps 



(2?. 



\K+1 



n + X+irn + 2K + 2 



(n + 2K + 2) 



I ^nK\i I ^^nK< n— 1 V uk' 



(59) 
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III. SOME SPECIAL CASES 



Let us study three special cases. At first, we study the s(s)-states {I — I — 0, i.e., n — ^l 
). It follows that the spin-orbit coupling term k,{k, + l)/r'^ — 0, and also the corresponding 
approximation to it in Eq. (27). In the presence of the exact spin symmetry limit (C^ = 0), 
the s-states {k — —1) energy equation becomes 



'^{2 + b){Mc'' + Er,,-,)Db n + 



2 



2{n + 5) 2 



(60) 
where 

+ (Mc^ + ^--0 > 1. (61) 



The upper- and lower-spinor wave functions are 

F„,_i(r) = M^^_,^^!l±^e-^'^'-{l - e-'^'-y ^F^ (-n, n + 2r^ + 25; 1 + 2r;; e""'') , (62) 

and 



na{n + 2r] + 26) ^-ar^s f ^-ar\v+i 

^■^'''-\Mc^ + E^,_,){l + 2rjy' e j 

X 2F1 (1 - n; n + 2 (77 + 5) + 1; 2 (1 + 77) ; e-"'") , (63) 

^ ^{2 + b){Mc' + E^,.^)Db n + S 
^ 2{n + S) 2 ' ^ ' 

where N'n-i is calculated in the Appendix C. Overmore, for the s-states {k — 1) in the exact 

pseudospin symmetry {Cps = 0), the energy equation (48) becomes 



with 



2 



2{n + 82) 



(65) 



with 



^^^ly-^\l'^~ ^'^''^ ~ ^"^^ j ' ^ = 0,l,2,---. (66) 

The wave functions given by Eqs. (52) and (54) become 
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= ■A/"ni ^^\|^^" e-''^"'-(l - e-"'-)^^ (-n, n + 2% + 2<52; 1 + 2%; e""^) , (67) 

with 

(2 + ^)^2 2n+ 1 + ^1 + 46^2 

% = , 7 , (68) 

2n + 1 + ^1 + 46^2 4 ' 

and 

+-^-(Mc^-E„.)(l + 2,,)(^-^ ) ) 

X 2F1 (l-n;n + 2(r;2 + ^2) + l;2(l + r;2);e-"'-), (69) 

where E^i 7^ Mc^. 

Second, we study the nonrelativistic hmit. In applying the following appropriate mapping 

— 7> 1^ to Eqs. (27) and (32), wc obtain the energy 
levels of the Schrodinger equation for any arbitrary orbital quantum number / as 



and the radial wave functions as 

^niir) = K /^'^^^^" e-^"-(l - e-°0' 2^1 {-n, n + 2?j + 26;l + 2?j; e""') , (71) 



with 



l + y(l + 2,r + ^i.6M>l,,-^^-iW -, (72) 



where = mi 7712/ (mi + m2) is the reduced mass of the two atoms and Afni is calculated in 
the Appendix C. 

Third, the case a ^ 0, the nonrelativistic limit of the bound state solutions of Eqs. (41) 
and (42) in exact symmetry limit {Cs — 0) can be obtained as 

£„, = Z,-|f , ^ V. P3) 

^ \l + 2n+^{l + 2lf + ^Dr1 

and the normalized wave functions: 
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where 

^?(lfzTTy L4(^,1^2i)V|z,.^-l|, (75) 

which are identical to the previous results obtained by the function analysis method [51], 
the NU method [52,53] and the EQR (cf. Eq. (29) in Ref. [54] obtained in i:)-dimensions) . 



IV. NUMERICAL RESULTS 

As in Ref. [56], we calculate the non-relativistic energy levels as function of various 
values of the parameter that controls the width of the potential well a = 0.05-0.30 fm~^ 
and equilibrium inter-nuclear distance = 0.4, 0.8 fm for various states with quantum 
numbers n and I. The atomic units h — M — 1 are used and the dissociation energy is set 
to D — 15 fmr^. In Table 3, we display our results with those ones calculated by using the 
conventional approximation scheme suggested by Greene and Aldrich [12] to deal with the 
centrifugal term /(/ + l)/r^ together with the values obtained from the numerical integration 
procedures based on the MATHEMATICA package programmed by Lucha and Schoberl [57]. 
Obviously, our results are closely approaching the ones obtained in [57] for both short range 
(small a) as well as for long range (large a) potential (see e.g., our recent works [39,58]). 
This means that our new approximated calculations using the approximation scheme (15) 
proposed recently by us provides much better approximation to the centrifugal term than 
that in Ref. [55] even for large a values (see, e.g., [40] and the references therein). At small 
values of a (Kratzer potential), we have also calculated the energy levels as a function of 
Te = 0.1 — 1.5 fm for various quantum numbers n and /. Hence, our numerical values of these 
energy levels are shown in Table 4. Overmore, Table 5 presents some numerical values for 
the eigeenergies of the Dirac valence states obtained from Eq. (27) with parameters M — 1.0 
fm~^, D — 15 fm~^, a — 0.1,0.3 fm"^ and Tg = 0.4,0.8 fm (exact spin symmetry case, 
i.e., = fm~^). We noticed that there are only positive energy bound state solutions in 
the spin symmetry limit. One can also see from Table 5 that there are degeneracies between 
the eigenstates (npi/2, np3/2), (^24/2, ^^(^5/2), {nf5/2,nfr/2), {ng7/2,ngg/2), etc. In fact, each 
of these eigenstates form a spin doublet. For instance, for any specific value of n, where 
n = 0, 1, 2, • • • , npi/2 with k = 1 is considered as the partner of nps/2 with k — —2. Thus, 
states that have the same radial n and orbital angular momentum I quantum numbers with 
j — 1 + 1/2 and j — I — 1/2 are degenerate [59]. 
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On the other hand, in Table 6, we present some numerical values for the eigeenergies of 
the Dirac hole states obtained from Eq. (50) with the previous choice of potential parameters 
for the case of pseudospin symmetry limit {Cps = 0, 5.0, —5.0 and —10.0 fm~^). One can 
see from Table 6 that there arc degeneracies between the eigenstates {nsii2i{n — 1)^3/2), 
{n - 1) /5/2), (^4/2, (n - 1) 5-7/2), (n/7/2, (n - 1) /19/2), etc. In fact, each of these 
eigenstates form a pseudospin doublet. For instance, for specific value of n = 1, l,si/2 with 
K — —1 is considered as the partner of 0^3/2 with k — 2. Thus, states that have pseudo 
orbital angular momentum I quantum numbers, radial n and n — 1 with j — I — 1/2 and 
j — I + 1/2, respectively, are degenerate [59]. 



V. CONCLUSIONS 



We have studied the approximate bound state solutions of the Dirac equation for the 
GMP model with any arbitrary spin-orbit k state under the conditions of the spin (pseu- 
dospin) symmetry limitation by means of the NU method including a new improved ap- 
proximation scheme to approximate the centrifugal (pseudo-centrifugal) barrier term. By 
setting E(r) (A(r)) to the spherically symmetric GMP, we have derived the solutions of the 
Dirac equation for the relativistic energy eigenvalues and associated two-component spinor 
wave functions for arbitrary spin-orbit k, state that provides an approximate solution to the 
spin- and pseudo-spin symmetry. The resulting solutions of the wave functions are being 
expressed in terms of the Jacobi polynomials (or hypergeometric functions) . We have shown 
that the present spin symmetry can be easily reduced to the non-relativistic solution once 
we set V{r) = S{r) (i.e., A(r) = or = 0). The non-relativistic hmits of our solution are 
obtained by imposing appropriate transformations and recalling k,{k, + 1) — )> /(/ + 1) in the 
spin symmetry limits. Furthermore, when a ^ 0, our results can be reduced to the well- 
known bound state solutions for the Kratzer potential model. If we choose the spin-orbit 
quantum number k — —1 {k — 1) ior spin (pseudospin) symmetry, the problem reduces to 
the exact s(s)-wave Dirac solution. We must point out that the numerical calculations for 
eigenenergy of the Dirac states involved in Eqs. (27) and (50) are sensitive to the choice 
of the parameters Cg, Cps, a, r^., D and M. It is found that the spin (pseudospin) limit 
Dirac eigenenergy valence (hole) states along with the two-component spinors are identical 
with the results obtained previously by other methods and works [15]. It is noticed that 
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these results are obtained in a much simpler fashion than Ref. [15]. Finally, Eqs. (27) and 
(50) can be used to evaluate the binding energies of the GMP for diatomic molecules such 
as CH, CO and N2 [2,56] in the relativistic framework with spin and pseudospin symmetry 
cases for any range of potential well a. Equations. (41) and (57) can be also used to evalu- 
ate the binding energy of the Kratzer potential in the relativistic framework with spin and 
pseudospin symmetry cases at small range potential well [a 0). 
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Appendix A: Generalization of the NU method 

We briefly review the Nikiforov-Uvarov essential polynomials, root, eigenvalues and wave 
functions (see Eqs. (5) - (11) of Ref. [40]) being expressed in terms of the parameters q 
(i = 1, 2, • • • , 13) together with A, B and C: 

(i) The relevant constants: 



C4 = 7^ 



(1 - Ci) , C5 = - (C2 - 2C3) , Cg = C5 A, 



C7 = 2C4C5 - S, C8 = C4 C, Cg = C3 (C7 C3C8) Ce, 

2 

Clo = ci 2c4 2i/c^ - 1 > -1, cii = 1 - ci - 2c4 H sfco, > -1, 

C3 

C12 = C4 + ^/c^ > 0, Ci3 = -C4 + — {^/c^ - C5) > 0. 



(Al) 



(ii) The essential polynomials: 



Tt{z) = C4 + C5Z - [{^/c^ + C3^/c^) Z - ^/c^ , 



(A2) 



k 



(C7 + 2C3C8) - 2y/c^, 



(A3) 



t{z) = 1 - (c2 - 2C5) Z-2 [(^ + C3^) Z - ^/C^ , 



(A4) 




2c3 - 2 ( -f C3 V^) < 0. 



(A5) 
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(iii) The energy equation: 

(C2 - C3) n + cgn^ - (2n + 1) C5 + (2n + 1) + 03^) + C7 + 2C3C8 + 2^5^ = 0. (A6) 

(iv) The wave functions: 

p{z)=z''''>{l-C3zy'\ (A7) 

= ^^12(1 - C3^)^l^ C12 > 0, Ci3 > 0, (A8) 

yn{z) = p(f^o,crO(^i _ 2c3z), c,o > -1, Cn > -1, (A9) 

Fn.{z) = AfnnZ'''{l - CsZ^' P^''^ - 2CsZ), (AlO) 

where Pn^''^\x), /x>— l,z/>— 1 and x G [—1,1] are the Jacobi polynomials with 

Pi"'^)(l - 2s) = i^lik (-n, l + a + P + n;a + l;s), (All) 

and is a normalization constants. Also, the above wave functions can be expressed in 
terms of the hypergeometric function as 

Fnn{z) = Nnnz''^ (1 - c^zf'^ ^F^ (-n, 1 + Cio + Cn + Ti' Cio + 1; c^z) , (A12) 

where C12 > 0, C13 > and z e [0, 1/C3] . 



Appendix B: Solving Quartic Energy Equation 

In order to solve the quartic equation (41b), the first step in the solution is to define the 

following variables 

3 11 3 11 

u = a2- -al; v = ai + -al - -0302; w = oq - — Og + ^^^302 - -0301; (Bl) 

which enable us to write down the related auxiliary cubic equation of the form: 

aEl. + bEl. + cEn. + d^O, (B2) 

with coefficients 

a = 1, 6 = c^ — (u^-4w), - — . (B3) 
2 ' 16 ^ ; ' 64 V y 

The next step is solving the cubic equation (B2) by defining the variables 

. c 16^ 2 6^ ifoc , 1 2 1 . 
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and then applying one of the following three cases: 

(i) When all three roots are real {h < 0): we define the variables 



. = yj^; k = cos- {-i),F = j=^H^ cos (^) , G = V^sin , 



(B5) 



where all the arguments in the trigonometric functions are in radians, to obtain the three 
possible roots of (B2) in the form 

and hence the four energies of the original quartic equation (41b) take the forms 



' :r(2)^(3) 4a4 ^ ^ „ Oiwiy 4 04 



V "'^ V nK / (2)W3) 404' V n« V n« / (2)W3) 4 04 

(B7) 

It is worth noting that whenever we have three real root we always choose the two non-zero 
(2) (3) 

roots; say, and E^^ in (B6) of the cubic equation, 
(ii) When only one root is real (/i > 0): the definitions 

R - --g + Vh] S=\^R] T = --g - Vh] U = \^] (B8) 

enable us to write down the three roots of (B2): 

and hence the four energies of Eq. (41b) are 



{2)W3) 4a4 /W2)^(3) 4 04 



(BIO) 

It is worth noting that whenever we have one real root and two complex roots we always 
choose the two complex roots. 
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(iii) When all three roots are real and equal {f — g — h — 0), then the roots of (B2): 



UK UK UK 



(Bll) 



and hence the four energies of the original quartic equation (41b) are 

■• = £7(3) ^ fi;(4) 
a 4a4' "'^ 4a4' "Va 4a4" 



£;(!) ^ -24'/? - - — • = ^ ^(4) ^ 2i'/^ ^ 



(B12) 



Appendix C: Normalization constants 

The normalization constant, Mni can be determined in closed form. We start by using 
the relation between the hypergeometric function and the Jacobi polynomials (see formula 
(8.962.1) in [49]): 

, ^, r(n + i/ + l) 

+ r(. + i) ' "=i> 

to rewrite the wave functions in (32) as 

Fn.{r) = Mn. J!:^^^';^^!:\ e-^-'^^{l - e-'^^Y^Pt-''^-'\l - 2e— ). (C2) 

From the normalization condition [un,i{'r)]'^ dr — 1 and under the coordinate change 
X — 1 — 2e~"^, the normalization constant in (B2) is given by 



n\T{2en, + 1) 
r(n + 2£„^ + l) 



' 'l — x\ /l + x\ fl+X 



[Pt"-''^-'\x)]' dx. 

(C3) 



The calculation of this integral can be done by writting 

1 + X ^ f 1 — X 



2 V 2 

and using the following two integrals (see formula (7.391.5) in [49]): 



/I 
(1 - x)"-' (1 + x)'^ [Pi''^''\x)Y dx ^ T 



,^^ r(n + ^+l)r(n + /. + l) 
n\ur(n + u + + ' 



(C4) 
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which is vahd for Re(i/) > and Re(//) > —1 and (see formula (7.391.1) in [49]): 

(1 - (1 + [pi'^'^' (-)] ' dx = 2-+^+^ (C5) 

which is vahd for Re(i/) > —1, Re(/x) > —1. Thus, the normahzation constant: 



1 



r(2£„« + 1) 



Q;e„«(n + + 5i) r(n + 2£„«; + l)r(n + 2£„« + 25i) 



2(n + (5i) 



n!r(n + 25i) 



1/2 



(C6) 



where < n, k < oo. In the s-wave {k, = —1) case, the above result is written explicitly as 



r(277 + i) 



ar]{n + r] + S)r{n + 2r] + l)r(n + 2r] + 2S) 



2{n + 5) 



nW {n + 25) 



1/2 



(C7) 



Also, the non-relativistic normalization constant is therefore obtained as 



r(2^+i) 



ar]{n + rj + S)T{n + 2t] + l)T{n + 2r/ + 26) 



1 1/2 



2{n + 5) 



n\T n + 26 



(C8) 



24 



[1] Z.H. Deng and Y.P. Fan, Shandong Univ. J. 7, 162 (1957). 

[2] L. Infeld and T.E. Hull, Rev. Mod. Phys. 23, 21 (1951); M.F. Manning and N. Rosen, Pliys. 
Rev. 44, 953 (1933). 

[3] J.W. Dabrowska, A. Khare and U.P. Sukhatma, J. Phys. A: Math. Gen. 21, L195 (1988); F. 
Cooper, A. Khare and U. Sukhatma, Phys. Rep. 251, 267 (1995); X. Zou, L.-Z. Yi and C.-S. 
Jia, Phys. Lett. A 346, 54 (2005). 

[4] A.D. S. Mesa and C. Quesne and Y.F. Smirnov, J. Phys. A: Math. Gen. 31, 321 (1998). 

[5] S. Fliigge, P. Walger and A. Weiguny, J. Molec. Spectroscopy 23 (3), 243 (1967). 

[6] C.S. Jia, Y. Li, Y. Sun, J.Y. Liu and L.T. Sun, Phys. Lett. A 311, 115 (2003). 

[7] C.S. Jia, X.L. Zeng and L.T. Sun, Phys. Lett. A 294, 185 (2002). 

[8] S. Codriansky, P. Cordero and S. Salamo, J. Phys. A: Math. Gen. 32, 6287 (1999). 

[9] Z. Rong, H.G. Kjacrgaard and M.L. Sage, Mol. Phys. 101, 2285 (2003). 
[10] S.H. Dong and X.Y. Gu, J. Phys.: Conf. Ser. 96, 012109 (2008). 
[11] R.L. Greene and C. Aldrich, Phys. Rev. A 14, 2363 (1976). 
[12] Z. Rong and M.L. Sage, Interdiscip. Sci. Comput. Life Sci. 1, 163 (2009). 
[13] S.L. Nyeo and LC. Yang, Phys. Rev. E 63, 046109 (2001). 

[14] C. Berkdemir, Nucl. Phys. A 770, 32 (2006); W.-C. Qiang, R.-S. Zhou and Y. Gao, J. Phys. 

A: Math. Theort. 40, 1677 (2007); A.D. Alhaidari, Phys. Rev. Lett. 89, 068901 (2002). 
[15] L.H. Zhang, X.P. Li and C.S. Jia, Phys. Scr. 80, 035003 (2009). 

[16] S.M. Ikhdair, C. Berkdemir and R. Sever, Spin and pseudospin symmetry along with orbital 

dependency of the Dirac-Hulthen problem, to appear in the Appl. Math. Comp. (2011). 
[17] S.M. Ikhdair and R. Sever, Appl. Math. Comp. 216, 911 (2010). 

[18] S.M. Ikhdair and R. Sever, Cent. Eur. J. Phys. 8 (4), 652 (2010); S.M. Ikhdair and R. Sever, 

Int. J. Mod. Phys. A 25 (20), 3941 (2010). 
[19] C. Berkdemir, A. Berkdemir and R. Sever, J. Phys. A: Math. Gen. 39, 13455 (2006). 
[20] A. Soylu, O. Bayrak and I. Boztosun, J. Phys. A: Math. Theor. 41, 065308 (2008); W.-C. 

Qiang, J.-Y. Wu and S.-H. Dong, Phys. Scr. 79, 065011 (2009); X.Y. Liu, G.F. Wei and C. Y. 

Long, Int. J. Theor. Phys. 48, 463 (2009). 
[21] A. Arda, R. Sever, C. Tezcan and H. Akgay, Chin. Phys. Lett. 27, 040306 (2010); O. Bayrak 



25 



and I. Boztosun, J. Phys. A: Math. Theort. 40, 11119 (2007). 
[22] C.-S. Jia, P. Guo, Y.-F. Diao, L.-Z. Yi and X.-J. Xie, Eur. Phys. J. A 34, 41 (2007); Y. Xu, 
S. He and C.-S. Jia, J. Phys. A: Math. Theor. 41, 255302 (2008); H. Akcay, J. Phys. A: Math. 

Theor. 42, 198001 (2009). 



[23 
[24; 
[25 
[26 
[27; 
[28 
[29 
[30 
[31 
[32 
[33 
[34; 
[35; 
[36 
[37; 

[38; 

[39 
[40 

[41 

[42; 
[43; 
[44; 
[45; 

[46; 



T. Chen, J.-Y. Liu and C.-S. Jia, Phys. Scr. 79, 055002 (2009). 

C.-S. Jia, J.-Y. Liu, L. He and L.-T. Sun, Phys. Scr. 75, 388 (2007). 

S.M. Ikhdair, J. Math. Phys. 51, 023525 (2010). 

O. Aydogdu and R. Sever, Phys. Scr. 80, 015001 (2009). 

C. Berkdemir and R. Sever, J. Phys. A: Math. Theor. 41, 045302 (2008). 

C.S. Jia, J.Y. Liu, L. He and L.T. Sun, Phys. Scr. 75, 388 (2007). 

J.N. Ginocchio, Phys. Rev. C 69, 034318 (2004). 

J.N. Ginocchio, Phys. Rev. Lett. 78, 436 (1997). 

J.N. Ginocchio, Phys. Rep. 414, 165 (2005); J.N. Ginocchio, Phys. Rep. 315, 231 (1999). 
J.S. BeU and H. Ruegg, Nucl. Phys. B 98, 151 (1975). 

P.R. Page, T. Goldman and J.N. Ginocchio, Phys. Rev. Lett. 86, 204 (2001). 

A. Arima, M. Harvey and K. Shimizu, Phys. Lett. B 30, 517 (1969). 

K.T. Hecht and A. Adler, Nucl. Phys. A 137, 129 (1969). 

A.D. Alhaidari, H. Bahlouli and A. Al-Hasan, Phys. Lett. A 349, 87 (2006). 

J. Meng, K. Sugawara-Tanabe, S. Yamaji, P. Ring and A. Arima, Phys. Rev. C 58, R628 

(1998). 

A.F. Nikiforov and V.B. Uvarov, Special Functions of Mathematical Physics, Birkhauser: 
Basel, 1988. 

S.M. Ikhdair, Phys. Scr. 83, 015010 (2011). 

S.M. Ikhdair, Eur. Phys. J. A 39, 307 (2009); S.M. Ikhdair and R. Sever, Phys. Scr. 79, 035002 
(2009). 

S.M. Ikhdair, Chem. Phys. 361, 9 (2009). 

W. Greiner, Relativistic Quantum Mechanics (Springer, Verlag, 1981). 
A.D. Alhaidari, H. Bahlouh and A. Al-Hasan, Phys. Lett. A 349, 87 (2006). 
C.-S. Jia, T. Chen, L.-G. Cui, Phys. Lett. A 373, 1621 (2009). 
Z.-Y. Chen, M. Li , C.-S. Jia, Mod. Phys. Lett. A 24 (23), 1863 (2009). 
Y.-F. Diao, L.-Z. Yi, T. Chen, C.-S. Jia, Mod. Phys. Lett. B 23, 2269 (2009). 

26 



[47] T. Chen, J.-Y. Liu, C.-S. Jia, Phys. Scr. 79, 055002 (2009), T. Chen, Y.-F. Diao, C.-S. Jia, 

Phys. Scr. 79, 065014 (2009).. 
[48] C.-S. Jia, Y.-F. Diao, L.-Z. Yi and T. Chen, Int. J. Mod. Phys. A 24 (24), 4519 (2009). 
[49] M. Abramowitz and LA. Stegun, Handbook of Mathematical Functions (Dover Publication, 

New York, 1972). 

[50] LS. Gradshtein and LM. Ryzhik, Tables and integrals, series and products (New York, Aca- 
demic, 1969). 

[51] S.M. Ikhdair and R. Sever, J. Mol. Struct.: Theochem 855, 13 (2008). 
[52] S.M. Ikhdair, Chin. J. Phys. 46, 291 (2008). 

[53] S.M. Ikhdair, Eur. Phys. J. A 40, 143 (2009); S.M. Ikhdair and R. Sever, Appl. Math. Comp. 
216, 545 (2010). 

[54] S.M. Ikhdair and R. Sever, J. Math. Chem. 45, 1137 (2009). 

[55] C. Berkdemir and Y.-F. Cheng, Phys. Scr. 79, 035003 (2009). 

[56] S.H. Dong and X.Y. Gu, J. Phys.: Conference Series 96, 012109 (2008). 

[57] W. Lucha and F.F. Schoberl, Int. J. Mod. Phys. C 10, 607 (1999). 

[58] S.M. Ikhdair and J. Abu-Hasna, Phys. Scr. 83, 025002 (2011). 

[59] R. Lisboa, M. Malheiro, A.S. de Castro, P. Alberto and M. Fiolhais, Phys. Rev. C 69, 024319 
(2004). 



27 



FIG. 1: The GMP model with D = 15 fm'^ for (a) various potential ranges a = 0.05, 0.15, 0.30 
fm~^ along with = 0.40 fm, and (b) various equilibrium inter-nuclear distances = 0.40, 0.80, 
1.20 fm along with a = 0.40 fm~^. 



FIG. 2: The upper and lower spinor wave functions, in the exact spin symmetry, for (a) n = 
(Opi/2; 0^3/2) spin doublet eigenstates with k = 1 and k = —2 and (b) ra = 1 (lpi/2) IP3/2) spin 
doublet eigenstates with k = 1 and k = —2. 



TABLE I: Specific values of the constants in the solution of the GMP under spin symmetry. 



Constant Constant 



ci = l 






C2 = 


1 




C3 = l 






C4 = 







C5 = -I 






C6 = 


1 
4 


1 + 4 ((2 + 6)^1 -^„.) 




C8 = 


— R 


C9 = i 


(1 + 2^^ + 46^1 




ClO = 


= 2?7i = 2iy^ 


cii = V^(1 + 2ac)' + 46^i 




C12 = 


= Vi 


Cl3 = 5 


(^l + y/(l + 2Acf + 46eij 


A = 


(2 + 6) - 


B = -2 


(^n« - ^1) - + 1) 


C = 


^riK 
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TABLE II: Specific values of the constants in the solution of the GMP under pseudospin symmetry. 



Constant 



Constant 



ci = l 

C3 = l 

c^ =2(EnK 
eg = 



C2 = l 
C4 = 



C6 



1 + 4 



((2 + ^)^2- 



k(k — 1) 



-Eri 



(1-2^)^ + 46^2 
cii = y^(1-2k)' + 46^2 

C13 = I fl ^ 



C8 = V2 

cio = 2r/2 = 2i\l En 

Cl2 = ??2 



(l-2Acr + 46e2j A={2 + h)^^-En 
B = -2(^En^-^^^ -k{k-1) C = ril = -Enn 
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TABLE III: The Schrodinger bound state energy levels Eni (in au) of the GMP as functions of a 
and Te for various states with D = 15 where fi = ji = 1. 

re = 0.4 Te = 0.8 

States a Present DG [53] LS [54] Present DC [53] LS [54] 



2p 


0, 


.05 


7.86080 


7.8606 


7.8628 


4.14088 


4.14068 


4.14208 
























0, 


.10 


7.95329 


7.95247 


7.95537 


4.21917 


4.21835 


4.2204 
























0, 


.15 


8.04508 


8.04322 


8.04724 


4.29737 


4.29552 


4.2987 
























0, 


.20 


8.13616 


8.13287 


8.13842 


4.37551 


4.37221 


4.3769 








re 


= 0.4 










re = 0.8 




0, 


.25 


8.22656 


8.22142 


8.22892 


4.45360 


4.44845 


4.4551 


States 


a 




Present DG [53] LS [54] 


Present 




0, 


.30 


8.31630 


8.30889 


8.31874 


4.53166 


4.52425 


4.5332 


4/ 


0. 


,05 


1 1 

J. ± . 


8908 


1 1 

-L ± . 


81 9^ 


1 1 

-L ± . 


.8209 


7.43469 


3p 


0, 


.05 


10.9978 


10.9976 


10.9998 


7.53279 


7.53258 


7.5350 




0. 


,10 


1 T 

± ± . 




1 1 

J. ± , 


. tJC/KJ 


1 1 

J. ± , 


.9981 


7.58636 




0, 


.10 


11.1626 


11.1617 


11.1647 


7.72475 


7.72393 


7.7271 




0. 


,15 


1 9 




1 9 




1 9 


.1718 


7.73559 




0, 


.15 


11.3242 


11.3224 


11.32647 


7.91516 


7.9133 


7.9177 




0. 


,20 








^991 


1 9 


.3421 


7.88251 




0, 


.20 


11.4828 


11.4795 


11.48513 


8.10400 


8.10071 


8.1066 


5p 


0. 


,10 


1 3 


.5421 


1 3 


5413 


1 3 

-LO. 


.5434 


11.3021 




0, 


.25 


11.6383 


11.6331 


11.64068 


8.29129 


8.28615 


8.2941 




0. 


,20 


1 3 


998Q 


1 3 


9957 


1 3 

-LO. 


.9301 


11.9132 




0, 


.30 


11.7907 


11.7833 


11.67565 


8.47703 


8.46962 


8.4799 


5d 


0. 


,10 


lO. 


.ouuo 


1 


.ou^o 


lO. 


.3075 


10.5201 


3d 


0, 


.05 


10.2160 


10.2154 


10.21651 


5.73974 


5.73913 


5.7404 




0. 


,20 


-LO. 




-LO. 




-LO. 


.6931 


11.0692 




0. 


.10 


10.3535 


10.351 


10.35409 


5.84574 


5.84327 


5.8465 


5/ 


0. 


,10 


1 3 


1475 


1 3 


.1426 


1 3 

-LO. 


.1478 


9.7966 




0. 


.15 


10.4893 


10.4837 


10.48992 


5.95061 


5.94505 


5.9515 




0. 


,20 


13, 


.5332 


13, 


.5134 


13, 


.5333 


10.2728 




0, 


.20 


10.6233 


10.6135 


10.62403 


6.05441 


6.04453 


6.0553 


5g 


0. 


,10 


13, 


.0379 


13, 


.0296 


13, 


.0379 


9.15212 




0, 


,25 


10.7557 


10.7403 


10.75645 


6.15720 


6.14177 


6.1582 




0. 


,20 


13, 


.4267 


13, 


.3938 


13, 


.42667 9.55246 




0, 


.30 


10.8864 


10.8642 


10.88719 


6.25904 


6.23682 


6.2601 


6p 


0. 


,10 


14, 


.0521 


14, 


.0513 


14, 


.0530 


12.2798 


Ap 


0. 


.05 


12.4976 


12.4974 


12.4992 


9.61301 


9.6128 


9.6156 


6d 


0. 


,10 


13, 


.9070 


13, 


.9045 


13, 


.9075 


11.7364 




0. 


.10 


12.6968 


12.696 


12.69851 


9.88351 


9.88269 


9.8862 


6/ 


0. 


,10 


13, 


.8111 


13, 


.8062 


13, 


.8113 


11.2448 




0, 


.15 


12.8883 


12.8865 


12.8901 


10.1485 


10.1467 


10.1514 


6g 


0. 


,10 


13, 


.7465 


13, 


.7383 


13, 


.7466 


10.8152 




0, 


.20 


13.0722 


13.0689 


13.0740 


10.4080 


10.4047 


10.4111 
























0. 


.25 


13.2484 


13.2433 


13.2501 


10.6619 


10.6568 


10.665 
























0. 


.05 


12.0983 


12.0977 


12.0989 


8.49334 


8.49272 


8.4948 
























0. 


.10 


12.2850 


12.2825 


12.2857 


8.70708 


8.70461 


8.7087 
























0, 


.15 


12.4664 


12.4608 


12.46715 


8.91774 


8.91218 


8.9194 
























0, 


.20 


12.6424 


12.6326 


12.64324 


9.12538 


9.11551 


9.1272 
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TABLE IV: The Schrodinger bound state energy levels E^i (in au) of the Kratzer potential as a 
function of for various states with D = 15 where h = /j, = 1. 



State/r 


e= 0.1 


0.4 


0.8 


1.0 


1.5 


2p 


1 S 

J- (J. 




7.76759 


4.06249 


3.21339 


2, 


0774Q 


3p 


14. 




10.8298 


7.33925 


6.26836 


4 


56776 


3d 


14. 


5192 


10.0766 


5.63256 


4.41694 


2, 


.74032 


4p 


14. 


,7319 


12.2908 


9.33707 


8.27299 


6, 


.40188 


4d 


14. 


7269 


11.9062 


8.27643 


7.04418 


5, 


.05045 


4/ 


14. 


724fi 


11.6401 


7.28048 


5.81633 




61468 


5p 


14. 


,8268 


13.0996 


10.6443 


9.65890 


7, 


.79162 


5d 


14. 


,8242 


12.8774 


9.94062 


8.80199 


6, 


,76420 


5/ 


14. 


,8230 


12.7278 


9.29975 


7.96875 


5, 


.69315 


59 


14. 


,8224 


12.6242 


8.73653 


7.18597 


4, 


.60125 


6p 


14. 


,8789 


13.5937 


11.5460 


10.6568 


8, 


.86976 


6d 


14. 


,8774 


13.4539 


11.0555 


10.0356 


8, 


.07050 


6/ 


14. 


,8768 


13.3616 


10.6191 


9.44444 


7, 


.25037 


6g 


14. 


,8764 


13.2985 


10.2430 


8.89964 


6, 


.42694 
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TABLE V: The eigeenergies of the Dirac valence states in units of fm for several values of n 
and K with the parameters M = 1.0 fm~^ and D = lb fm~^ in the case of exact spin symmetry 
limit {Cs = fm~^). We have set ^ = c = 1. 



I 




' *'? 


K <0 


Tllj A 1 i A In 
J=l+1/ 2 




n 




> nLj^i_ 






— 10 fm^l 

\J. l.\J J 1 1 V ^ 


' e 


— 40 fm 

\J.'-t\J J lib 














1 




0, 


-2 




5.5791076 


0, 


1 




5.5791076 






n 


— o 


U"5/2 


6.8118605 


0, 


9 


U"3/2 


D.oiloOUO 


3 







—4 




8.0171073 


0, 


3 


Oft;/n 
'^J5/2 


8 0171073 


4 




0, 


-5 


OOq /o 


9.1025175 


0, 


4 


OQ'7 /o 


9.1025175 


1 




1 


_2 


-'■^'3/2 


8.1823677 


1, 


1 


-'-/^l/2 


8 1823677 


2 




1, 


-3 


Idt; /o 


8.8815340 


1, 


2 


liiq /o 


8.8815340 


3 




1, 


-4 


J 7/2 


9.6603105 


1, 


3 


J 5/2 


9.6603105 


4 




1, 


-5 


lOo /o 


10.4200196 1, 


4 




10.4200196 


OL 


— fi fm~^ 

— U.OU Jill ^ 




— n 40 fm 














1 







-2 




5.7078594 


0, 


1 




5 7078594 


9 
z 




n 


Q 
— O 


U«5/2 


6.9646771 


0, 


9 
z 


U«3/2 


fi QR/LR771 


3 







—4 


Of-T/O 


8.2121326 


0, 


3 


0/"c;/n 

'^J5/2 


8 2121326 


4 




0, 


-5 


0(7q/o 


9.3506414 


0, 


4 


0(77/9 


9.3506414 


1 




1 

-■-5 


_2 


-'-K3/2 


8.4626850 


1, 


1 


1 Til In 

^P\/2 


8 4626850 


2 




1, 


-3 


Idc; /o 


9.1762544 


1, 


2 


Idq /o 


9.1762544 


3 




1, 


-4 


1 r? 


9.9831712 


1, 


3 


lf^/9 


9.9831712 


4 




1, 


-5 


IOq /o 


10.7812870 1, 


4 


l0'7 /O 


10.7812870 


Cv 
\X 


— 10 fm~^ 


' e 


— 80 fm 














1 







-2 


ODq /o 


3.6831690 


0, 


1 


0?9i /o 


3.6831690 







n 


O 


U«5/2 


4.3378367 


0, 


9 


0*^3/2 




3 




0, 


-4 


0/7/2 


5.0775317 


0, 


3 


0/5/2 


5.0775317 


4 




0, 


-5 


05y/2 


5.8291637 


0, 


4 


O57/2 


5.8291637 


1 




1, 


-2 


lP3/2 


5.8388616 


1, 


1 


lPl/2 


5.8388616 


2 




1, 


-3 


14/2 


6.2180443 


1, 


2 


14/2 


6.2180443 


3 




1, 


-4 


1/7/2 


6.6999906 


1, 


3 


1/5/2 


6.6999906 


4 




1, 


-5 


lfy/2 


7.2334098 


1, 


4 


I57/2 


7.2334098 
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TABLE VI: The eigenenergies of the Dirac hole states in units of fm for several values of n and 
K with the parameters M = 1.0 fm~^ and D = 15 fm~^ in the case of the pseudospin limit. We 
have set ^ = c = 1. 



/ 




n, 


K < 


/J -X 


En,K<0 


En,K<0 


n 


- 1, 


K>0 {l + 2,j + 1) En-l,n>0 




a 


= 0.10 /m"\ 




= 0.40 /m 




Cps = 


Cps = 5.0 








C^pg 


Cps = 5.0 


1 




1, 


-1 




7.1975980 


9.0681299 


0, 


2 




7.1975980 


9.0681299 


2 




1, 


-2 


lP3/2 


7.9184956 


9.7118773 


0, 


3 


0/5/2 


7.9184956 


9.7118773 


3 




1, 


—3 


l«5/2 


8.6773995 


10.3429382 


0, 


4 




8.6773995 


10.3429382 


4 




1, 


—4 


J- 77/2 


9.3887565 


10.9088876 


0, 


5 


Ohq/o 


9.3887565 


10.9088876 


1 




2, 


— 1 


2si/2 


8.7054710 


10.1707692 


1, 


2 


Of z 


8.7054710 


10.1707692 


2 




2, 


—2 


2p3/2 


9.1768973 


10.6114449 


1, 


3 


1/5/2 


9.1768973 


10.6114449 


3 




2, 


—3 


^"5/2 


9.7101038 


11.0728353 


1, 


4 


1^/7/9 


9.7101038 


11.0728353 


4 




2, 


—4 


2/7/2 


10.2352574 11.5024429 


1, 


5 


Iha/o 

a 1 £1 


10.2352574 11.5024429 


a 


= 0.30 /m~\ 


J^e 


= 0.40 /m 




Cps — 


Cps — —5.0 








Cps — 


Cps — —5.0 


1 




1, 


-1 


lSl/2 


7.4924218 


6.2513272 


0, 


2 


Odrt /o 


7.4924218 


6.2513272 


2 




1, 


-2 


lP3/2 


8.2327097 


6.9536060 


0, 


3 


Of5/9 


8.2327097 


6.9536060 


3 




1, 


-3 


lc^5/2 


9.0263198 


7.7574142 


0, 


4 


C7 / / Z 


9.0263198 


7.7574142 


4 




1, 


—4 


1/7/2 


9.7816289 


8.5599024 


0) 


5 


O/iq/9 


9.7816289 


8.5599024 


1 




2, 


— 1 


2si/2 


9.1118917 


8.0615654 


1, 


2 




9.1118917 


8.0615654 


2 




2, 


—2 


2^3/2 


9.5965089 


8.5254559 


1, 


3 


If 5/2 


9.5965089 


8.5254559 


3 




2, 


-3 


2^5/2 


10.1539358 9.0879925 


1, 


4 


1^/7/9 


10.1539358 9.0879925 


4 




2, 


—4 


2/7/2 


10.7110039 9.6762585 


1, 


5 


y/z 


10.7110039 9.6762585 


a 


= 0.10 /m~\ 


J^e 


= 0.80 fm 




Cps — 


Cps = -10.0 








Cps — 


Cps = -10.0 


1 




1, 


-1 


lSl/2 


4.9972498 


2.8369925 


0, 


2 


o/z 


4.9972498 


2.8369925 


2 




1, 


-2 


lP3/2 


5.4281366 


3.1294486 


0, 


3 


0/5/2 


5.4281366 


3.1294486 


3 




1, 


-3 


lc^5/2 


5.9435202 


3.5245459 


0, 


4 


057/2 


5.9435202 


3.5245459 


4 




1, 


-4 


1/7/2 


6.4876355 


3.9879463 


0, 


5 


0/i9/2 


6.4876355 


3.9879463 


1 




2, 


-1 


2Sl/2 


6.3778620 


4.3337279 


1, 


2 


1^^3/2 


6.3778620 


4.3337279 


2 




2, 


-2 


2^'3/2 


6.6705985 


4.5479115 


1, 


3 


1/5/2 


6.6705985 


4.5479115 


3 




2, 


-3 


2^5/2 


7.0449439 


4.8443318 


1, 


4 


I57/2 


7.0449439 


4.8443318 


4 




2, 


-4 


2/7/2 


7.4613765 


5.2011464 


1, 


5 


1^9/2 


7.4613765 


5.2011464 
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